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ABSTRACT 

An  atteapt  is  ns  da  to  show  theoretically  that  tha  drag  caused  by  a 
turbulent  boundary  layer  on  a  snoot h  flat  plate  vibrating  with  low  amplitude 
transverse  standing  naves  la  the  sane  as  for  a  non-vibrating  plate  provided 
the  wave  celerity  of  the  conponent  traveling  naves  la  greater  than  about  3.5 
tines  the  free  strean  velocity.  Experiments  conducted  in  a  water  tunnel  with 
a  rectangular  cross-section  hsving  one  vibrating  wall  shewed  no  Measurable 
change  in  boundary  layer  velocities, as  compared  to  non-vibrating  conditions, 
for  frequencies  ranging  fron  15  to  122  cps  and  frse  strean  velocities  rang¬ 
ing  fron  2.U  to  15  fpe.  Tha  lowest  wave  oalerity  to  free  strean  velocity 
ratio  tested  was  2.5U* 
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TURBULENT  BOUNDARY  LAYER  FLOW 
OVER  A  FLAT  PLATE  VIBRATING 

WITH  TRANSVERSE  STANDING  WAVES 

I.  INTRODUCTION 

Fluid  drag  on  a  flat  plate  in  steady  motion  is  fairly  easy  to  esti¬ 
mate.  Recently,  increased  attention  has  been  given  to  flow  past  wavy  bound¬ 
aries.  This  interest  was  first  directed  to  the  problem  of  the  generation  of 
surface  waves  by  wind.  Benjamin  [1]  summarises  the  work  of  previous  re¬ 
searchers  in  this  field,  notably  that  of  Miles  [2]. 

It  is  generally  established  that  drag  is  inposedon  the  flow  by  the 
wave  by  a  sheltering  effect  wherein  the  pressure  on  the  windward  side  of  the 
wave  is  higher  than  on  the  leeward  side.  This  effect  can  be  explained  con¬ 
veniently  by  the  phenomenon  of  separation.  Miles,  however,  by  taking  viscos¬ 
ity  into  account,  showed  that  sheltering  can  occur  without  separation.  Ben¬ 
jamin  in  his  paper  generalizes  the  work  of  Miles  and  computes  a  sheltering 
effect  for  a  laminar  boundary  layer  profile.  He  proposes  that  the  method  is 
extendable,  with  limitations,  to  turbulent  boundary  layer  profiles.  He  re¬ 
ports  agreement  with  the  work  of  Motzfeld  [3]  who  experimentally  discovered 
a  sheltering  pressure  distribution  without  separation  for  flew  over  fixed 
waves  with  steepnesses  of  0.025. 

More  recently,  tests  by  Kramer  [U]  indicated  that  stabilization  of 
laminar  boundary  layers  could  be  accomplished  by  applying  to  the  surface  a 
compliant  dissipative  coating  in  which  traveling  waves  are  generated  by  the 
flow.  A  discussion  and  analysis  of  this  phenomenon  is  also  given  by  Benjamin 
[5].  In  the  above  situations  it  is  important  to  realise  that  the  wave  celer¬ 
ity  c  is  always  less  than  the  main  stream  flow  velocity. 

Another  problem  of  waves  on  boundary  surfaces  involves  the  plates 
of  ship  hulls.  Vibrations  within  the  ship  cause  standing  waves  to  be  excited 
in  the  plates.  The  objective  of  the  research  on  which  this  report  is  based 
was  to  determine  if  these  standing  waves  influence  the  boundary  layer  to  the 
extent  of  causing  an  additional  drag.  Tests  were  conducted  in  a  horizontal 
water  tunnel  having  a  rectangular  cross-section  8  in.  deep  and  12  in.  wide. 
A  flexible  steel  plate  12  in.  wide  and  U5  in,  long  composed  part  of  the  upper 

* 

Numbers  in  square  brackets  refer  to  the  List  of  References  on  p.  27. 
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surface.  This  plate  was  vibrated  with  frequencies  ranging  from  15  to  122  cps. 
Velocity  traverses  upstream  and  downstream  of  the  vibrating  section  were  made 
for  several  frequencies  within  this  range  and  for  core  velocities  ranging  from 
2.U  to  15  fps.  In  these  tests  the  boundary  layer  appeared  completely  unaf¬ 
fected  by  the  state  of  vibration.  Analysis  based  on  the  work  by  Benjamin  [1] 
on  traveling  waves  has  been  adapted  to  standing  waves  by  superposition  to 
confirm  this  result. 

This  research  was  carried  out  under  the  Bureau  of  Ships  Fundamental 
Hydromechanics  Research  program,  SR-009-01-01 , administered  by  the  David  Tay¬ 
lor  Model  Basin.  Grateful  acknowledgment  is  expressed  to  the  members  of  the 
Laboratory  staff  who  participated  in  this  program.  Particular  acknowledgment 
is  expressed  to  K.  Yalamanchili  who  aided  in  the  preliminary  stages  of  the 
project,  to  Z.  S.  Tarapore  who  collected  much  of  the  data, and  to  Carol  Takyi 
who  prepared  the  manuscript  under  the  direction  of  Loyal  Johnson.  This  study 
has  been  conducted  under  the  general  direction  of  Dr.  Lorenz  G.  Straub,  Di¬ 
rector  of  the  Laboratory. 

II.  THEORETICAL  DESCRIPTION  OF  THE  FLOW 

A.  General  Theory  for  Traveling  Waves 

In  his  paper,  Benjamin  [1]  analysed  the  action  of  shearing  flews 
passing  over  a  boundary  surface  disturbed  by  a  continuous  train  of  harmonic 
waves.  The  waves  travel  along  the  flexible  but  stationary  surface  with  a 
celerity  c.  His  analysis  parallels  that  used  in  the  study  of  the  stability 
of  laminar  boundary  layers.  In  the  stability  study  a  simplification  of  the 
Navier-Stokes  equation,  known  as  the  Orr-Summerfeld  equation,  is  used  to  an¬ 
alyse  the  behavior  of  small  amplitude  traveling  wave  disturbances  assumed  to 
exist  wit  hint  he  boundary  layer  flow.  The  relationship  between  wave  celerity, 
wave  length,  and  amplitude  attenuation  is  determined  for  different  Reynolds 
numbers.  Conditions  favoring  a  negative  attenuation  are  assumed  to  be  un¬ 
stable. 

The  Orr-Sumamrfeld  equation  also  arises  in  this  analysis  of  flow 
over  a  traveling  surface  disturbance.  However,  in  this  case  the  wave  length, 
amplitude ,  and  celerity  are  given  a  priori  as  a  property  of  tbs  boundary.  The 
character  of  the  disturbance  to  the  flow  caused  by  the  wave  is  determined  to 
find  out  if  the  wall  stresses  are  such  as  to  cause  an  extra  drag  on  the  sur¬ 
face  .  The  extra  drag  is  caused  by  a  sheltering  effect  in  which  the  pressure 
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distribution  is  out  of  phase  with  the  wave  form  so  that  the  pressures  on  the 
downstream  side  of  the  wave  are  less  than  those  on  the  upstream  side.  This 
is  not  to  be  eonfusedwith  the  amplifying  effect  caused  by  pressure  distribu¬ 
tions  in  phase  with  the  wave  that  is  important  in  stability  studies.  The 
salient  points  of  the  analysis  are  presented  in  the  remaining  paragraphs  of 
this  subsection. 

Two  transformations  are  introduced  at  the  beginning  of  the  analysis 
to  bring  the  boundary  conditions  into  a  more  tractable  form.  First,  to  a- 
chieve  a  steady  flow  situation  the  co-ordinate  system  is  referenced  to  the 
wave  celerity.  The  wave  form  is  thereby  fixed  in  space  but  the  surface  it¬ 
self  and  the  flow  are  given  an  additional  velocity  component  -c.  Second, 
curvilinear  orthogonal  co-ordinates  are  used,  corresponding  to  the  stream 
function  and  the  potential  function  for  potential  flow  over  the  wavy  surface. 
The  co-ordinate  17  is  directed  across  the  flow  and  the  co-ordinate  £  is 
directed  along  the  flow.  The  £,17  system  is  correlated  to  a  rectangular 
x,  y  system  by  the  complex  transformation. 


K  *  i>? 


x  +  iy  -  iae 


ik(x  ♦  iy) 

f 


(1) 


where  X  is  the  wave  length  and  a  is  the  wave  amplitude .  The  equation  for 
the  boundary  surface  is  simply  ij  ■  0.  With  the  flow  transformed  onto  the 
£,17  plane  the  boundary  is  a  flat  surface. 

The  stream  function  for  the  disturbed  flow  ^  (£,17  )  is  assumed  to 
have  two  components.  The  first  ^(17)  is  the  stream  function  for  the  main 
shear  flow  assuming  no  wave,  i.e.  17  -  y.  It  is  given  by 


where  11(17)  is  the  undisturbed  velocity.  The  second  oomponent  ^(£,17  )  is 
the  periodic  harmonic  perturbation  caused  by  the  wave  and  is  written  for  later 
convenience  in  the  form 

+p  ■  (F(t)  *  W*)  -  c).-*1*]  (3) 


where  F(ij)  is  a  residue  function  to  be  computed.  The  complex  notation  takes 
into  account  a  possible  phase  shift  in  reference  to  the  wave  form.  The 
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remaining  analysis  requires  that  the  wave  steepness  be  small  so  that  the  terms 
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k  and  ak  are  of  second  order  magnitude. 


The  velocity  of  the  disturbed  flow  to  first  order  in  ak  oan  be 
described  by  components  in  the  5  and  rj  directions  using  the  respective 
equations. 


11(5,77  )  •  -  U  -  c  ♦  a[f'  ♦  u'e“k7>]  e1*5 

)  — J*  -  -ika[F  ♦  (0  -  c)e"kT?]  e^5 
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where  J  is  the  Jacobian  of  the  co-ordinate  transformation.  Here  again  the 
complex  notation  allows  for  a  phase  shift  even  though  only  the  real  parts  of 
u  and  v  are  to  have  meaning. 

The  two-dimensional  Havier-Stokes  equations  for  laminar  flow,  when 
written  in  terms  of  the  stream  function,  yield  upon  elimination  of  the  pres¬ 
sure  terms  the  following  expression  for  F. 

(U  -  c)(F*  -  k2F)  -  u"f  -  -i  [F1*  -  2k2?2  *  kVje’1^  (5) 

.mi  . 

assuming  parallel  main  flow  (U  and  U  »  0) .  This  is  the  Orr-Summerfeld 
equation  which  is  basic  to  the  theory  of  boundary  layer  stability.  (See 
Schlicting  [6],  page  316.)  For  reasonably  high  Reynolds  numbers,  the  viscous 
terms  on  the  right  are  small  and  negligible  except  in  special  friction  layers. 
One  of  these  occurs  at  the  boundary  and  another  occurs  in  the  region,  if  any, 
where  U(i?)  approaches  c.  In  our  problem,  c  will  be  shewn  to  be  larger 
than  U  for  all  practical  cases  so  that  this  second  layer  will  not  exist. 

The  function  Tty)  can  be  resolved  handily  into  two  components 
which  are  additive.  The  first  t±ty)  is  an  inviscid  function  which  satis¬ 
fies  (5)  when  the  right  hand  terms  are  neglected  to  give 

(0  -  c)(^  -  (6) 


The  second  i8  *  viscous  function  which  is  negligible  every¬ 

where  except  near  the  boundary  and  which  satisfies  the  equation 

w  -  •>*;  -  -1  -*•  4T 


(7) 
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This  equation  ia  obtained  from  Eq.  (5)  by  neglecting  the  terns  containing 
U  and  k^.  The  terms  containing  these  values  are  small  at  the  boundary. 

Both  ^  and  ^  *ust  approach  sero  as  The  other  bound- 

ary  conditions  are 

^(o)  ♦  ■  e  and  ^(o)  ♦  j(^(o)  “  -®'(o)  (8) 


Solutions  of  ^  and  ^  approach  sero  at  infinity  will  each  contain 
one  undetermined  coefficient.  The  coefficient  for  ^  can  be  obtained  frcai 
the  combined  boundary  condition 


^(o)  -  -U*(o) 


(9) 


The  coefficient  for  ^  follows  from  and  Eq.  (8), 

The  pressure  on  the  boundary  as  related  to  the  pressure  at  infin¬ 
ity  can  be  obtained  by  integrating  the  simplified  Navier-Stokes  equation  along 
a  5  line  from  17-0  to  17  -  cO .  The  function  ^  is  significantly  large 
for  17  up  to  about  -y-  but  the  region  of  significant  i8  ▼•ry  small* 
The  function  ^  ia  therefore  negligible  inthis  integration  and  the  pressure 
oan  be  expressed  by 

- (D 

po  ’  pc0  *  Po#lk4»  Po  "  P^2  /  (D  "  c)*ldr>  ^°) 

J0 

Shear  due  to  the  wave  action,  on  the  other  hand,  is  almost  entirely  attrib¬ 
utable  to  jfg  and  is  given  by 


To  "  ^2  (°)ielk*  (11) 

P0  end  rQ  are  seen  to  be  harmonic  with  respect  to  £  but  possibly  out  of 
phase  with  the  wave  form. 

The  extra  drag  induced  by  the  wave  action  is  computable  from  these 
boundary  stresses.  First  from  Iq.  (10)  it  follows  that  the  real  and  Imaginary 
components  of  ^  will  give  respectively  real  and  imaginary  components  of 
Fq.  The  real  component  P0p  will  give  a  negative  pressure  distribution 
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phased  with  the  wave  fora.  This  distribution  will  contribute  zero  net  fora 
drag.  The  imaginary  component  P0^  gives  a  pressure  distribution  -jj-  out 
of  phase  with  the  wave  and  a  net  drag  does  result.  This  is  the  so-called 
sheltering  effect. 

The  resulting  drag  coefficient  is  given  by  the  expression 


akP. 


„  _  n 

cw  "TT 


p»: 


oO 


(12) 


The  harmonic  shear  distribution,  of  course,  contributes  no  net  drag. 

The  function  ^  would  ordinarily  be  real  since  Eq.  (6)  has  only 
real  parts.  However,  the  function  ^  complex  and  through  the  combined 
boundary  conditions  of  Eq.  (9)  a  complex  character  can  be  transferred  to 
Therefore,  although  ^  is  essentially  inviscid  in  origin,  its  final  form 
is  influenced  by  the  viscid  function  The  effect  of  ^  on  1®  *B 

though  a  "displacement  thickness"  were  created  along  the  boundary.  The  thick¬ 
ness  varies  harmonically,  but  out  of  phase,  with  the  wave  form  producing  a 
falsely-shaped  boundary  for  the  inviscid  flew. 

In  effect  more  important  in  boundary  layer  stability  and  wave  gen¬ 
eration  is  the  complex  character  that  ^  receives  due  to  the  singularity 
of  the  point  where  U  »  c.  As  mentioned  before,  however,  this  situation  is 
not  important  to  this  paper. 


B.  A  Traveling  Harmonic  Wave  with  No  Main  Flow 

The  simplest  model  to  which  this  theory  can  apply  is  the  one  in 
which  there  is  no  main  shear  flow.  This  case  is  not  treated  ty  Benjamin  [1]. 
However,  investigation  of  this  problem  will  bring  out  the  character  of  the 
analysis,  especially  as  it  applies  to  the  present  study.  U(iy)  now  is  ident¬ 
ically  sero  and  Eqs.  (6),  (7),  and  (8)  have  the  exact  solutions 


(13) 

(HO 
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The  expression  for  ^  is  complicated,  but  if  the  variable  kv/c 
is  snail  £q.  (lU)  reduces  to 

^(-l)  ■  o.-k1  (IS) 

Since  both  k  and  v  are  snail,  this  equation  must  hold  except  for  the  very 
slowest  wave  celerities  (c  <  l).  The  velocities  associated  with  obtain* 

able  from  Eq.  (U) ,  are 


u(S,»7 )  ■  -akce"*^  cos  k(5  -  ct) 

,  N  _kn  (16) 

vU,i})  ■  +akce  "  sin  k(£  -  ct) 

In  these  equations  the  wave  form  celerity  has  been  restored.  This  is  exact* 
ly  the  potential  flew  under  a  moving  wave.  The  wall  pressure  distribution 
from  Eq.  (10)  has  the  form 

Pc  -  pco  *  -p(kc)2  -j-  cos  k(£  -  ct)  (17) 

Figure  1  shows  the  potential  flew  for  this  problem  and  the  resulting  veloc¬ 
ities  and  stress  distributions.  The  depth  at  which  the  potential  velocity  is 
1  per  cent  of  the  boundary  velocity  is  0.73X . 

The  viscid  solution  ^(ij)  is  also  recognisable.  It  is  equivalent 
to  the  boundary  layer  on  a  flat  oscillating  plate  (see  Schlictlng  [6],  page 
67  ).  The  u  component  of  velocity  becomes 

)  -  akc  axp  (  ij)  coe(  *  k(*  -  ct))  (18) 

The  wall  shear  due  to  the  viscid  solution  is  given  by 

T0  ■  -pa(kc) ®in  (k*  ♦  -J-  *  ket)  (19) 

The  velocity  pattern  and  shear  for  this  function  are  shown  as  well  in  Fig.  1. 
The  thickness  over  which  the  viscid  solution  is  significant  has  been  oalled 
the  wall  friction  layer.  Its  thickness  may  be  estimated  by  the  magnitude  of 
the  quantity  ^2v/kc  which  is  for  most  conditions  vexy  small. 
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The  above  analysis  shows  that  the  flow  pattern  beyond  the  wall  fric¬ 
tion  layer  is  very  accurately  represented  by  the  usual  potential  flow  except 
for  eases  where  kv/c  approaches  1  or  larger.  This  usually  requires  very 
low  wave  celerities.  The  viscosity  is  shown  to  be  directly  influential  only 
in  a  very  narrow  friction  layer  next  to  the  boundary. 

The  possibility  of  a  sheltering  effect  due  to  a  zone  of  separation 
in  the  lee  of  the  waveform  is  definitely  not  present  in  this  case.  A  fluid- 
filled  cavity  following  the  wave  would  have  to  be  transported  over  the  sta¬ 
tionary  boundary  at  a  velocity  c.  This  situation  is,  obviously,  impossible. 

C.  A  Standing  Harmonic  Wave  with  No  Hain  Flow 

The  second  model  worthy  of  consideration  is  the  case  of  a  standing 
harmonic  wave  as  shown  in  Fig.  2.  A  standing  wave  with  amplitude  a,  wave 
length  X  ,  and  frequency  f  can  be  produced  by  the  superposition  of  two 
traveling  wave  trains  with  equal  and  opposite  celerities.  These  waves  must 
each  have  an  amplitude  of  £a,  a  wave  length  of  X  and  celerities  c  equal 
to  plus  and  minus  Xf.  It  is  useful  to  note  here  that  kc  *  2xf.  The  flow 
pattern  for  a  standing  wave  can  be  obtained  by  the  corresponding  superposi¬ 
tion  of  the  traveling  wave  flow  patterns.  This  is  entirely  acceptable  within 
the  limitations  of  small  wave  steepnesses  and  the  neglect  of  second  order 
terms. 

Figure  2  shows  tbs  result  of  superimposing  the  flow  patterns.  When 
the  plate  is  in  the  null  position  the  potential  streamlines  and  velocity  pat¬ 
terns  are  exactly  the  same  as  for  traveling  waves.  In  the  loop  position  the 
inviscid  velocities  are  everywhere  zero.  It  is  interesting  to  note  that  the 
viscous  velocities  reinforce  st  the  nodes  but  cancel  at  the  loops.  Although 
it  is  not  shown  in  the  figure,  the  viscous  velocities  do  not  all  become  zero 
when  the  plate  is  in  the  loop  position  due  to  a  time  phase  shift. 

The  pressure  distribution  hss  a  standing  wave  pattern  in  phase  time- 
wise  and  anti-phase  spacewise  with  the  wave  form.  In  the  case  of  the  travel¬ 
ing  wave,  a  sheltering  effect  could  only  be  achieved  if  the  quantity  kv/c 
approached  1.  In  the  present  case  the  sheltering  effect,  if  any,  is  auto¬ 
matically  cancelled  out .  This  result  would  be  suggested  also  by  the  symmetry 
of  the  flow  pattern. 
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There  is,  of  course,  no  separation  possible  in  this  case  either. 
The  particles  near  the  boundary  are  confined  to  small  orbital  paths  and  do 
not  accumulate  in  a  zone  of  separation. 


D.  A  Traveling  Wave  with  a  Turbulent  Boundary  Layer 


Turbulent  boundary  layers  would  seem  to  be  applicable  to  the  gen* 
eral  theory  as  long  as  the  friction  layer  lies  within  the  laminar  sublayer, 
or  at  least  not  too  far  beyond  it. 

Most  of  the  analysis  in  the  preceding  paragraphs  cones  from  Ben¬ 
jamin  [1].  For  the  case  of  boundary  layer  flow  solutions  of  Eqs.  (6),  (7), 
and  (8)  are  not  obtained  exactly.  Simplification  of  Eq.  (7)  can  be  accom- 
plished  by  assuming  the  velocity  of  the  boundaiy  layer  U  to  be  linear  as 

i 

in  the  laminar  sublayer,  so  that  U  « 77  U  (0).  All  the  coefficients  in  Eq. 
(7)  are  thereby  constant  with  rj  .  Equation  (7)  in  this  new  form  has  been 
solved  numerically  in  connection  with  boundary  layer  stability  studies.  Tab¬ 
ular  values  of  the  function  ^  are  presented  for  instance  in  Holstein  [7]. 
The  variable  77  ,  however,  is  replaced  by  a  new  variable,  z,  defined  by 

me  /  cU  (o)  \ 

z  -  mn  -  — | - ,  where  m  ■  ( - )  (20) 

0  (o)  \  v  / 

Notice  that  on  the  boundary  zQ  ■  -mc/D '  ( o) .  The  function  jfg  decreases 
very  rapidly  with  increasing  z  or  if  so  that  the  wall  friction  layer  is 
again  very  thin.  It  is  found  to  be  negligible  by  the  time  i)  reaches  l/m. 

Of  special  value  in  Eq.  (9)  is  the  Tietjens  function  D(-zq)  de¬ 
fined  by 


D(-z0) 


(21) 


The  real  and  imaginary  parts  of  0  are  tabulated  in  Holstein  [7],  and  graphed 
in  Benjamin  [1].  Both  components  of  D  approach  zero  as  -*o  becomes  fair¬ 
ly  large.  In  ary  event,  it  is  0(1)  or  less  for  practical  values  of  -a  . 

o 

A  simple  approximate  inviscid  solution  to  Eq.  (6)  is 
^  -  +B(c  -  U)e"kl? 


(22) 
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It  holds  fairly  well  as  long  as  no  singularity  occurs  as, for  instance,  where 
U  «  c.  The  coefficient  B  comes  from  consideration  of  Eq.  (9)  and  has  the 
form 


B 


1  ♦  D(-»0)/»0 
1  ♦  D(-so)/zo  +  D(-s0)/m 


(23) 


If  the  term  -zq  is  0(D(-zQ)),  D  in  Eq.  (23)  will  dominate  and  ^  will 
be  ccaiplex ,  giving  rise  to  a  sheltering  effect.  However,  if  -zq  is  much 
larger  than  D,  B  will  approach  the  real  value  1,  ^  will  be  real,  and 
there  will  be  no  sheltering  effect.  Inspection  of  the  function  D  shows 
that  for  -zq  >  8,  |o(-zo)|  is  less  than  0.25  and  decreases  to  approach 
zero  as  -aQ  increases.  Therefore  for  -zQ>  8  sheltering  would  be  negli¬ 
gible. 

It  is  a  relatively  simple  matter  to  show  that  ~*0>  8  is  met  pro¬ 
viding  c  >  SUqq*  From  the  study  of  turbulent  boundary  layers  (see  Schlict- 
ing  [6],  page  U32.)  the  velocity  gradient  at  the  wall  can  be  related  to  the 
core  velocity  by  the  equation 


<j'(°)  -  /2v  (2U) 

i 

where  cf  is  the  local  drag  coefficient  dependent  on  the  Reynolds  number. 
Introducing  Eq.  (21;)  into  the  expression  defining  zq  and  eliminating  m 
gives  the  relationship 


(25) 


According  to  Schlicting,  for  Reynolds  number  greater  than  10^  the  quantity 

i  i 

c.  is  less  than  0.0055*  Substitution  of  -s  "8,  v  ■  10  ,  and  c_  ■ 

^  O  a 

0.0055  into  Eq.  (21;)  gives  c  ■  U.l#^.  It  follows  therefore  that  for  tur¬ 
bulent  boundary  layers  the  general  theory  indicates  a  negligible  sheltering 
effect  provided  the  wave  celerity  is  greater  than  about  five  times  the  main 
stream  velocity. 


It  remains  to  determine  a  criterion  which  insures  that  the  general 
theory  Itself  is  applicable  to  a  particular  situation.  The  restriction  is 
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that  the  wall  friction  layer,  which  extends  to  about  17  ■  ,  should  lie 
within  the  region  of  linear  velocity.  According  to  Benjamin  this  region  ex- 
tends  to  77  -  10v/V  even  though  the  usual  expression  for  the  extent  of  thB 
laminar  sublayer  is  rj  ■  5v/Vq.  Equating  the  expressions  for  these  two  re¬ 
gions  and  introducing  again  Eq.  (2U)  gives  the  expression 


hOOOv 


Substituting  v  »  10“^  and  c'  -  0.005$  produces  the  result  that  c  -  5.2U-. 

1  Ou 

A  wave  celerity  greater  than  five  times  the  main  stream  velocity 
would  seem  to  be  an  adequate  criterion  to  insure  that  there  be  no  sheltering 
due  to  the  effects  of  the  wall  friction  layer.  Apparently  the  wave  length 
and  wave  amplitude  are  not  directly  significant.  The  effect  of  the  boundary 

layer  thickness  is  reflected  in  the  effect  of  Reynolds  number  on  c '  For 

9  *  * 

R  >  10  ,  c^  is  less  than  0.001U  and  the  criterion  becomes  c  greater  than 

2-6%- 

Sheltering  caused  by  a  zone  of  separation  following  the  wave  form 
is  also  Impossible  tinder  the  above  criterion.  The  reasoning  is  the  same  as 
given  in  the  case  of  a  traveling  wave  with  no  flow.  A  convenient  point  of 
view  from  which  to  observe  this  situation  is  with  the  reference  velocity  ad¬ 
justed  so  that  the  wave  form  is  stopped.  With  this  reference  it  can  be  seen 
that  there  is  no  area  of  flow  reversal  as  long  as  c  is  greater  than  U  _ . 

Tbs  only  remaining  mechanism  by  which  the  traveling  wave  could  in¬ 
fluence  the  main  flow  is  interference  in  the  turbulence  pattern.  This  pos¬ 
sibility  is  completely  unlikely  as  long  as  the  wall  friction  layer  lies  with¬ 
in  the  laminar  sublayer  as  required  by  the  above  criterion.  The  oscillating 
movements  of  the  potential  flow  alone  produce  no  strong  shears  that  could  in¬ 
fluence  the  turbulence. 

The  equation  for  the  wall  pressure  distribution  for  a  turbulent 
boundary  layer  may  be  obtained  by  substituting  Eq.  (22)  for  with  B  «  1 
into  Eq.  (10)  for  pQ  to  get 
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P0  _poO  “  “I*2*®1*5  fU-  -^^“)2c2e“k7?d7?  (27) 

JQ  0 

can  be  substituted  for  U (17)  without  materially  affecting  the  results, 
especially  if  c  is  large.  Upon  integration  of  Eq.  (27)  and  restoration  of 
the  wave  velocity  the  expression  for  pressure  is 

p0 "  P<D  "  • (1  -  -~-)2  p(kc)2  Hr 008  k(5  *  ct)  (28) 

which  is  comparable  to  Eq.  (17)  for  the  case  of  no  flow. 

E.  A  Standing  Wave  with  a  Turbulent  Boundary  Layer 

The  research  of  this  report  involved  a  case  of  turbulent  boundary 
layer  over  a  boundary  surface  vibrating  with  standing  waves.  A  criterion  for 
the  existenoe  of  a  sheltering  effect  for  this  case  can  be  obtained  by  super¬ 
imposing  the  appropriate  turbulent  boundary  layer  traveling  wave  situations. 
It  is  first  necessary,  however,  to  reconsider  the  previous  case  of  traveling 
waves  with  a  turbulent  boundary  layer  but  with  the  flow  velocity  and  the  wave 
celerity  opposite  in  direction  to  each  other.  The  variable  a  must  be  rede¬ 
fined  compared  with  Eq.  (20)  to  give 


s 


">T?  ’ 


mo 

U*(o) 


(29) 


where  the  definition  of  n  is  unchanged.  Since  a  and  U*(o)  in  this  case 
are  numerically  negative,  the  relationship  between  s  and  ^t*)  and  rj 
is  exactly  as  before.  In  order  to  make  use  of  existing  tables  and  graphs, 
the  TLetjens  function  should  be  redefined  as 

D(-t0)  -  ♦  m^2(o)/^(o)  (30) 


When  these  adjustments  are  made,  the  criterion  previously  established  can  be 
applied  equally  well  to  this  reversed  case. 

For  analysis,  the  turbulent  boundary  layer  passing  over  a  standing 
wave  must  be  separated  into  two  equivalent  boundary  layers,  each  moving  over 
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oppositely  directed  traveling  waves.  These  equivalent  boundary  layers  will 
each  contain  half  the  undisturbed  main  flow  so  that  all  velocities  are  half 
of  the  undivided  flow  velocities.  The  important  distances  from  the  wall  such 
as  the  thickness  of  the  boundary  layer  and  the  thickness  of  the  laminar  sub¬ 
layer  are  unchanged  by  this  arbitrary  division.  The  wall  velocity  gradient 
will  be  Ju  (o)  where  the  symbol  refers  to  the  undivided  flow.  A  boundary 
layerwith  a  free  stream  velocity  equal  to  ,  developing  naturally,  would, 
however,  have  quite  different  proportions.  The  thickness  of  the  laminar  sub¬ 
layer  would  be  roughly  twice  that  of  the  undivided  flow  and  the  wall  velocity 
gradient  would  be  only  roughly  £u  (o) .  The  naturally  developing  boundary 
layer  is  unimportant  to  this  analysis,  however,  since  the  criterion  to  be 
developed  depends  only  on  the  character  of  which  in  turn  depends  only 
on  the  conditions  at  the  boundary  and  not  on  the  method  of  development. 

To  find  a  relationship  between  c  and  1)^  for  the  undivided  flow, 
it  is  only  necessary  to  refer  to  Eqs.  (25)  and  (26).  In  Eq.  (25)  c  must 
be  adjusted  so  that  -z  remains  constant  at  8.  Therefore  c  must  vary  as 

i  o 

the  root  of  c_.  This  relationship  holds  for  Eq.  (26)  as  well.  Substituting 

f  * 

JU  (o)  and  _  into  Eq.  (2U)  shows  that  the  local  drag  coefficient  for 

QJ  I  I 

the  divided  flow  is  2cf .  Substituting  2cf  and  into  Eq.  (25)  or 

Eq.  (26)  shows  that  c  must  be  50^  or  about  3.50^.  Therefore  for 

Reynolds  numbers  based  on  distance  from  the  leading  edge  greater  than  about 

105  there  is  negligible  sheltering  due  to  the  wave  disturbance  provided  o 

9  1 

is  greater  than  3.50^ .  For  R  greater  than  10  tbs  lower  values  of 
require  that  c  need  only  be  greater  than  21)^ . 

Separation  as  a  source  of  sheltering  is  impossible  in  this  case  as 
well  for  the  reasons  given  earlier.  It  can  be  further  remarked  that  separa¬ 
tion  is  caused  by  a  sustained  adverse  pressure  gradient.  Pressure  gradients 
are  present  in  this  case,  of  course,  but  at  any  one  spot  they  alternate  rap¬ 
idly  between  favorable  and  adverse. 

The  wall  pressure  distribution  for  a  turbulent  boundary  layer  over 
a  standing  wave  is  given  by  the  appropriate  superposition  of  the  traveling 
wave  pressure  distributions  as  represented  by  Eq.  (28)  only  with  the  free 
stream  velocity  given  as  .  The  resulting  expression  is 


1U 


u 


-  -  [1  +  )2]p(2nf)2  -£L  sin  k5  Bin  kct 


po  “  pcO  "  '  lA  '  '  2c  7  2n 


+  »yP-  p(2nf)2  -4-^-  cos  kg  cos  kct 
ZC  tn 


(31) 


III.  CHARACTERISTICS  OF  A  VIBRATING  PLATE 

The  behavior  of  a  plate  vibrating  with  transverse  standing  waves 
can  be  analyzed  by  the  classical  theory  of  vibrations  in  elastic  beams.  The 
phenomenon  can  be  thought  of  as  characterized  by  the  interaction  of  the  in¬ 
ertia  forces  due  to  the  mass  of  the  plate  with  the  elastic  forces  due  to  bend¬ 
ing  of  the  plate.  When  the  plate  forms  the  boundary  of  a  fluid  the  inertial 
effects  of  the  fluid  must  also  be  considered.  In  the  case  of  liquids  the 
elasticity  of  the  fluid  is  ordinarily  negligible. 

When  an  infinitely  long  plate  of  unit  width  and  depth  d  is  vi¬ 
brating  with  simple  harmonic  standing  waves  of  amplitude  a,  wave  length  X* 
and  frequency  f  the  inertial  force  due  to  the  mass  of  the  beam  is  given  by 

■  Spda(2nf)2  sin  Cos  2 lift  (32) 

where  p  is  the  density  of  water  and  S  is  the  specific  gravity  of  the  ma¬ 
terial  of  the  beam.  The  distribution  of  wall  pressure  for  the  case  of  a 
standing  wave  with  no  flow  is  given  in  Fig,  2  to  be 

Pc  ■  -  P  -g—-  a(2nf)2  sin  cos  2 aft  (33) 

The  negative  of  the  pressure  nay  be  thought  of  as  the  inertial  force  of  the 
water  per  unit  of  plate  area.  The  negative  is  required  since  positive  pres¬ 
sures  exert  forces  opposite  in  direction  to  y.  Equation  (33)  holds  for  cases 
with  flow  provided  is  small  as  seen  in  Eq.  (31) .  Comparison  of  Eqs. 
(32)  and  (33)  shows  that  the  term  -g^-  corresponds  to  d.  The  effect  of  the 
fluid  is  therefore  to  augment  the  mass  of  the  plate  by  the  mass  of  a  layer  of 
water  of  depth  .  The  frequency  of  vibration  can  now  be  computed  using 
the  classical  theory  (see,  for  example,  Timoshenko  [8],  page  331.)  to  give 


f  - 


Oh) 


12(pSd  ♦  p\/2m) 
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where  E  ie  Young's  modulus  of  elasticity.  Equation  (3U)  can  be  put  in  more 
useful  form  using  the  relationship  c  ■  f  X  to  give,  in  terms  of  d/X,  the 
relationship 


c 


E 


12pS(l  +\  /2ndS) 


(35) 


Figure  3  shows  a  plotting  of  wave  celerity  versus  the  ratio  dA 
based  on  Eq.  (35) .  The  physical  constantswere  evaluated  for  steel  and  water. 
The  term  for  the  added  mass  of  the  water  is  included  in  the  lower  curve  but 
neglected  in  the  upper  curve .  This  chart  emphasises  the  strong  dependence  of 
wave  celerity  on  the  stiffness  of  the  plate. 

The  above  remarks  apply  to  plates  or  beams  infinitely  long.  A  plate 
of  finite  length,  however,  pin  supported  at  its  ends,  will  vibrate  as  an  in¬ 
finitely  long  plate  with  nodes  located  at  the  supports.  The  reason  is  that 
an  infinitely  long  plate  has  aero  bending  moment  and  maximum  shear  at  the 
nodes  as  should  occur  at  a  pin-supported  end.  A  given  plate  will  have  a  num¬ 
ber  of  modes  of  vibration.  The  fundamental  mode  with  only  one  loop  is  worse 
from  the  point  of  view  of  this  study  since  it  has  the  lowest  equivalent  wave 
celerity.  A  plate  of  length  L  vibrating  in  its  fundamental  mode  will  have 
a  wave  length  of  2L. 


Practical  aspects  of  design  naturally  limit  the  slenderness  of  ac¬ 
tual  plates  and  by  doing  so  limit  the  equivalent  wave  celerities  c  of  res¬ 
onant  standing  waves.  The  midspan  deflection  A  of  a  uniformly  loaded  beam, 
simply  supported  at  its  end,  is  given  by  the  expression 


L  E  38U  d 


(36) 


where  w  is  the  loading  per  unit  length.  For  example,  steel  plate  loaded 
with  one  atmosphere  of  pressure  will  have  a  ratio  of  about  0.01  for  an 
L/d  ratio  of  50.  The  corresponding  wave  celerity  is  175  fps.  For  a  L/d 
ratio  of  100  the  ratio  is  0.08  and  the  corresponding  wave  celerity  is 
72  fps.  The  latter  values  of  sag  seen  far  beyond  the  limits  of  acceptable 
design. 

The  steepness  of  standing  waves  is  also  limited  by  design  condi¬ 
tions.  Figure  3  shows  the  relationship  between  the  d/X  ratio  andthe  wave 
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steepness  ratio  a/X  as  governed  by  two  different  design  criteria.  The  first 
is  the  "extreme  fiber  stress"  in  the  deflected  plate  and  the  second  is  the 
fluctuation  in  vail  pressure.  The  vave  steepness  as  related  to  plate  thick¬ 
ness  is  shown,  by  way  of  example,  for  an  extreme  fiber  stress  of  20,000  psi, 
a  common  structural  value.  This  curve  shows  that  except  for  extremely  flex¬ 
ible  plates  the  wave  steepness  is  restricted  to  the  order  of  1  per  cent  or 
less.  The  other  curve  on  this  chart  shows  the  wave  steepness  that  will  cause 
a  pressure  fluctuation  of  plus  and  minus  one  atmosphere.  For  this  condition 
cavitation  vouldoccur  at  the  loops  if  the  reference  pressure  is  atmospheric. 
This  limitation  is  seen  to  be  very  severe  for  the  larger  d/X  ratios. 

IV.  EXPERIMENTAL  PROGRAM 

The  purpose  of  the  experiments  was  to  determine  whether  or  not 
standing  wave  vibrations  on  a  smooth  boundary  plate  would  disturb  a  turbulent 
boundary  layer  sufficiently  to  cause  an  extra  induced  drag.  It  was  recognised 
that  the  method  of  testing  must  be  such  that  an  effect  would  be  detected  even 
if  it  were  very  small.  It  was  necessary  that  the  measurements  be  as  funda¬ 
mental  and  direct  as  possible  in  order  that  the  turbulence  of  the  flow  and 
the  vibrations  in  the  tunnel  vouldnot  mask  the  overall  effect  that  was  sought. 
It  was  impossible  to  attach  any  instrument  to  the  plate  itself  as  the  vibra¬ 
tions  of  the  plate  would  make  the  readings  meaningless.  Wall  shear  measur¬ 
ing  devices  such  as  a  shear  plate  dynamometer  or  a  shear  tube  of  the  preston 
type  were  considered  for  location  Just  downstream  of  the  plate.  But  it  has 
been  reported  that  the  wall  shear  adjusts  itself  very  rapidly  to  the  boundary 
conditions  sothat  the  effect  of  the  vibrating  plate  upstream  would  not  likely 
be  sensed  by  such  devices. 

The  effect  of  vibration  on  the  boundary  layer  is  measured  most  di¬ 
rectly  and  completely  by  velocity  traverses  taken  upstream  and  downstream  of 
the  vibrating  section.  Comparison  of  these  velocity  profiles  for  the  dif¬ 
ferent  conditions  of  vibration  provides  an  immediate  test  for  the  presence 
of  any  effect.  Furthermore,  introduction  of  this  data  into  the  momentum  e- 
quation  would  provide  a  quantitative  estimate  of  the  change  in  drag.  It  is 
important,  of  course,  that  this  data  be  very  accurately  measured  and  that  all 
extraneous  effects  be  reduced  to  a  minimum. 
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V.  TEST  APPARATUS 

The  apparatus  on  which  the  experiments  were  performed  le  shown  dia- 
grammatically  in  Fig.  U.  Figure  5  contains  photographs  of  the  equipment  as 
well.  The  vibrating  plate  water  tunnel,  fabricated  for  this  research,  is  a 
modification  of  a  permanent  facility  of  the  Laboratory.  Prior  to  alteration, 
the  water  tunnel  consisted  of  a  12  in.  by  18  in.  rectangular  test  section 
drawing  water  from  a  vertical  stilling  tank  10  ft  high  and  U  ft  in  diameter. 
The  tank  is  supplied  from  the  Mississippi  River  by  a  12  in.  diameter  penstock 
with  a  control  valve  at  the  inlet  to  the  tank.  The  pressure  in  the  tunnel 
was  controlled  by  a  throttling  gate  at  the  discharge  end.  The  portions  of 
the  earlier  tunnel  used  in  the  present  facility  are  shown  in  bold  lines  in 
Fig.  lu 

It  is  essential  in  boundaiy  layer  studies  that  the  turbulence  of 
the  incoming  flow  be  as  low  as  possible.  The  stilling  tank  provided  was  ac¬ 
tually  too  small  for  its  purpose  considering  that  there  is  a  Uo-ft  pressure 
head  on  the  penstock  which  jets  water  into  the  tank  at  $0  fps.  The  cross- 
sectional  area  of  the  tunnel  was  changed  to  12  in.  by  8  in.  to  reduce  the 
discharge  required  and  make  the  tank  more  effective .  Screens  and  flow  straight- 
eners  were  provided  to  reduce  the  turbulence  further  to  a  seemingly  acceptable 
level. 

The  actual  tunnel,  constructed  of  lucite,  was  made  with  a  cross- 
section  12  in.  square.  The  upper  8  in.  of  depth  formed  the  tunnel  proper. 
'The  lower  L  in.  was  separated  by  a  slotted  wall  to  form  a  plenum  chamber  with 
its  own  pressure  control.  The  plenum  chamber  had  two  purposes.  One  was  to 
provide  some  control  over  the  longitudinal  pressure  gradient  in  the  tunnel 
and  the  other  was  to  increase  the  depth  of  the  tunnel  to  the  order  of  so 
that  fluctuating  pressures  and  velocities  would  be  negligible  at  the  lower 
wall. 

The  condition  of  zero  pressure  gradient  along  the  wall  requires,  of 
course,  a  constant  core  velocity.  The  effect  of  the  boundary  layer  is  to  re¬ 
duce  the  effective  cross-section  dimensions  of  the  tunnel  by  the  amount  of 
the  displacement  thickness,  causingths  core  velocity  to  increase.  Discharge 
must  therefore  be  bled  off  through  the  slotted  wall  as  the  boundaiy  layer 
displacement  thickness  increases.  The  slot  discharge  per  unit  of  wall  area 
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Is  proportional  to  the  free  stream  velocity  and  the  lengthwise  rate  of  in¬ 
crease  of  the  displacement  thickness.  The  rate  of  increase  of  the  displace¬ 
ment  thickness  is  approximately  proportional  to  the  distance  from  the  lead¬ 
ing  edge  to  the  power  of  -  1/5*  (See  Schlicting  [6],  page  U32.)  The  slot¬ 
ted  portion  started  21  in.  downstream  of  the  leading  edge  of  the  transition 
and  extended  for  5  ft  of  the  tunnel.  Thus  the  required  slot  discharge  was 
highest  at  the  upstream  end  and  tapered  off  slightly  downstream. 

The  slot  discharge  per  unit  area  of  wall  is  proportional  to  the 
slot  velocity  and  the  slot  porosity.  The  velocity  through  the  Blots  is  con¬ 
trolled  by  the  pressure  drop  across  them.  The  slot  velocity  at  different 
points  along  the  wall  is  influenced  by  lengthwise  variations  in  the  pressures 
in  the  main  tunnel  and  in  the  plenum  chamber.  The  pressure  variations  in 
the  plenum  were  minimised  by  making  the  chamber  sufficiently  large  so  that 
the  velocity  heads  associated  with  the  plenum  flowwere  everywhere  negligible* 
The  effect  of  pressure  variation  in  the  main  flow  can  be  minimized  if  the 
pressure  drop  across  the  wall  is  of  a  much  larger  order  of  magnitude  than  the 
variations.  A  good  criterion  is  that  the  velocity  through  the  slots  be  made 
equal  to  the  main  flow  velocity.  Then,  assuming  the  slots  as  orifices,  a  1 
per  cent  variation  in  tunnel  pressure  will  produce  only  0.5  per  cent  varia¬ 
tion  in  slot  velocity. 

With  the  slot  velocity  equal  to  the  main  stream  velocity  the  proper 
distribution  of  slot  discharge  is  accomplished  if  the  porosity  varies  direct¬ 
ly  with  the  rate  of  growth  of  the  displacement  thickness.  On  the  basis  of 
preliminary  calculations  a  slot  porosity  distribution  was  used  with  porosi¬ 
ties  ranging  from  0.0055  upstream  to  0.00U5  downstream. 

The  vibrating  portion  of  the  tunnel  was  a  stainless  steel  plate  12 
in.  wide,  U5  in.  long,  and  1/8  in.  thick,  composing  part  of  the  upper  wall. 
It  was  mounted  completely  independent  of  the  tunnel  and  connected  to  it  only 
by  a  very  thin  flexible  plastic  seal.  The  plate  was  supported  at  each  end 
by  a  heavy  structural  steel  frame  fastened  directly  to  the  basement  floor  of 
the  laboratory.  Special  care  was  required  to  align  the  ends  of  the  plate  to 
the  remainder  of  the  lucite  upper  wall  to  insure  that  there  would  not  be  a 
vertical  discontinuity  in  the  upper  boundary  and  that  there  would  be  no  phys¬ 
ical  contact  with  the  rest  of  the  tunnel.  The  end  supports  were  pin  connected 
so  that  the  plate  would  vibrate  naturally  as  an  infinite  plate  with  nodes  at 
the  ends.  The  nodes  are  points  of  zero  bending  moment  but  maximum  shear. 
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The  supporting  frame  had  to  be  rigid  enough  to  withstand  the  large  vertical 
harmonic  end  reactions  supplied  by  the  plate  and  to  give  a  minimum  of  verti¬ 
cal  displacement. 

A  1/8  in.  thick  plate,  US  in.  long,  has  little  resistance  to  internal 
pressure.  Two  intermediate  vertical  restraints  were  provided  to  keep  the 
plate  from  bulging  out  of  its  plane.  Their  position  was  adjustable  so  that 
they  could  be  located  at  intermediate  nodal  points.  Even  so,  the  internal 
pressure  had  to  be  carefully  controlled.  This  was  accomplished  by  a  head 
loss  gate  at  the  discharge  end. 

Various  types  of  electronic,  electro-mechanical,  and  mechanical 
vibrators  were  investigated  before  deciding  on  the  type  adopted  for  these 
tests.  The  exciter  used  was  a  pneumatic  vibrator  manufactured  by  the  Hartin 
Engineering  Company  of  Neponset,  Illinois  under  the  trade  name  "Vibrolator" . 
A  portion  of  one  is  shown  in  one  of  the  photos  of  Fig.  5.  It  consists  simply 
of  an  enclosed  circular  raceway  on  which  a  steel  ball  is  free  to  roll.  A  Jet 
of  high  velocity  air  strikes  the  ball  once  each  revolution  to  maintain  its 
circular  motion.  The  vibrator  transmits  a  two-directional  harmonic  vector 
force  to  any  object  on  which  it  is  bolted.  In  the  apparatus  the  vibrator  is 
mounted  to  the  vibrating  plate  at  mid-span  through  a  linkage  arrangement  that 
transmits  only  the  vertical  component  of  the  vector  force.  The  vibrator  is 
suspended  from  the  supporting  frame  by  lew  tension  springs  which  counterbal¬ 
ance  the  weight  of  the  vibrator  and  the  plate  but  transmit  little  of  the  vi¬ 
brational  force  to  the  frame. 

This  type  of  exciter  has  several  advantages  over  the  others  inves¬ 
tigated  for  this  program.  It  is  a  great  deal  more  compact  and  far  less  ex¬ 
pensive  than  electronic  exciters  and  has  the  further  advantage  that  the  vec¬ 
tor  force  increases  with  the  frequency  of  vibration  instead  of  decreasing  as 
is  the  usual  case.  In  addition  it  provides  a  pure  harmonic  force  with  a  fre¬ 
quency  that  is  controlled  simply  by  adjusting  the  air  pressure.  Its  main 
drawback  is  that  the  vector  force  is  dependent  entirely  cm  the  frequency. 
This  shortcoming  can  be  avoided  by  using  a  series  of  vibrators,  each  with  a 
different  vector  force-frequency  relationship. 

Only  two  models  in  the  Vibrolator  series  were  used  in  these  exper¬ 
iments.  One  was  the  model  CV35  which  had  a  0.383-lb  ball  and  a  raceway  pro¬ 
viding  a  0.625-in.  eccentricity  to  the  center  of  gravity.  The  vector  force 
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for  this  model  is  21*0  lb  at  100  cps.  The  other  model  was  the  CV19  with  a 
0.0622-lb  ball  and  a  0. 688-in.  eccentricity.  Its  vector  force  at  100  cps  is 
1*3  lb.  The  vector  force  varies  with  the  square  of  the  frequency. 

The  vibrational  displacements  of  the  plate  were  measured  with  e 
simple  displacement  meter  designed  at  the  Laboratory  employing  a  linear  var¬ 
iable  differential  transformer  produced  by  Schaevitz  Engineering  of  Camden, 
New  Jersey.  The  output  of  the  meter  was  fed  to  a  paper  strip  pen  recorder 
and  a  cathode  tube  oscilloscope.  The  pen  recorder  was  most  useful  for  meas¬ 
uring  frequencies  and  the  oscilloscope  for  measuring  displacement.  The  os¬ 
cilloscope  was  also  very  useful  for  observing  the  shape  of  the  wave  with  re¬ 
spect  to  time. 

A  series  of  flattened  stagnation  tubes  with  varying  tip  dimensions 
were  fabricated  in  order  to  measure  the  boundary  layer  velocities.  The  di¬ 
mensions  of  the  tubes  are  shown  in  Fig.  6.  All  but  one  of  the  tubes  were  made 
from  1/8- in.  diameter  stainless  steel  tubing  with  a  1/lS-in.  internal  diam¬ 
eter.  The  smallest  tube  of  tne  series,  however,  had  a  tip  made  of  l/L&~in. 
stainless  steel  tubing  with  a  l/32-in.  internal  diameter.  The  wall  thickness 
at  the  tip  was  reduced  to  0.002  in.  and  the  tip  was  flattened  to  a  thickness 
of  0.0085  in.  The  stems  of  all  the  tubes  of  the  series  were  made  from  lA-in. 
diameter  stainless  steel  tubing. 

Glands  were  provided  in  the  upper  wall  at  three  stations  for  in¬ 
sertion  of  the  stagnation  tubes.  The  gland  for  Station  1  was  2  in.  upstream 
of  the  plate,  for  Station  2  was  6  in.  downstream  of  the  plate,  and  for  Sta¬ 
tion  3  was  12  in.  downstream  of  the  plate.  Station  3  was  used  only  for  meas¬ 
uring  the  core  velocity.  The  actual  measurements  were,  of  course,  taken  at 
the  tipof  the  stagnation  tube,  which  is  2.5  in.  upstreamof  the  stem  portion. 

A  dial  gage  set  against  a  block  clasped  to  the  stagnation  tube  stem 
was  used  to  measure  the  distance  of  the  probe  from  the  wall.  In  addition,  a 
gold-plated  electrical  contact  was  set  flush  into  the  upper  lucite  wall  just 
over  the  stagnation  tube  tip.  An  electric  circuit  was  set  up  so  that  a  light 
turned  on  when  the  probe  just  touched  the  wall.  This  system  proved  very  sens¬ 
itive  and  provided  accurate  readings  for  the  traverse  distance. 

The  static  pressure  required  for  the  determination  of  the  velocity 
heads  was  measured  by  two  wall  taps,  each  located  on  the  upper  wall  two  inches 
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on  either  side  of  the  centerline  at  the  same  longitudinal  section  as  the  stag¬ 
nation  tube  tip.  Additional  vail  taps  were  Installed  along  the  back  vail  at 
mid-depth  In  the  8 -in.  channel  in  order  to  surrey  the  longitudinal  pressure 
profile. 


VI.  TEST  RESULTS 

The  measurements  describingthe  modes  of  vibration  of  the  plate  are 
shovn  on  Fig.  8.  The  plate  could  be  made  to  vibrate  vith  almost  a  continuous 
range  of  frequencies  by  adjusting  the  pressure  of  the  compressed  air  driving 
the  exciters.  However,  except  for  three  resonant  conditions  the  oscillations 
were  very  small  and  usually  negligible  except  under  the  exciter.  Careful  tun¬ 
ing  to  the  resonant  conditions  produced  the  very  good  vibration  patterns  shown 
in  Fig.  8.  Once  tuned,  the  vibrator  frequency  would  lock  in  phase  and  remain 
stable  even  for  fairly  substantial  changes  in  air  pressure.  Data  are  shown 
for  core  velocities  of  2.U  and  15  fps. 

Resonance  with  three  loops  and  vith  five  loops  was  obtained  vith 
the  larger  of  the  tvo  vibrators.  The  frequencies  for  the  two  conditions  are 
lli.8  cps  and  57  cps  and  the  amplitudes  under  the  exciter  are  0.050  and  0.011 
in.  Only  the  seven  loop  resonance  was  obtained  with  the  smaller  vibrator. 
The  frequency  and  amplitude  for  this  condition  are  122  cps  and  O.OOl*  in.  The 
amplitude  of  the  different  loops  along  the  plate  can  be  seen  to  vary.  This 
is  due  largely  to  the  concentrated  mass  of  the  exciter  mount  and  to  a  lesser 
degree  to  the  conditions  at  the  end  supports  where  some  movement  was  observed. 

The  wave  pattern  vith  respect  to  time  is  shovn  in  Fig.  8  on  photo¬ 
graphs  of  portions  of  the  pen  recorder  paper.  The  higher  frequencies  are  ap¬ 
proaching  the  limit  of  the  recorder  but  the  sinusoidal  character  of  the  wave 
is  clearly  evident .  The  pen  tracings  have  been  retouched  in  spots  to  improve 
the  reproduction  for  the  report.  The  frequencies  reported  were  measured  with 
this  recorder  but  the  amplitudes  were  obtained  from  oscilloscope  patterns. 

The  variation  in  amplitude  along  the  plate  caused  slight  variation 
in  wave  length,  but  the  nominal  values  are  30,  18,  and  12  6/7  in.  The  cor¬ 
responding  values  of  wave  celerity  c  are  36.9,  85.5.  and  133  ffcs  and  the 
d/x  values  for  the  l/8-in.  plate  are  0.001*2,  0.0071*,  and  0.0096.  The  wave 
celerity  data  are  plotted  on  the  chart  of  c  versus  d/X  on  Fig.  3.  The 
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points  are  seen  to  lie  below  the  theoretical  line.  The  difference  is  prob¬ 
ably  attributable  to  the  added  mass  of  the  exciter  mount.  Also  plotted  on 
Fig.  3  ere  the  ware  steepness  values,  0.0016,  0.00061,  and  0.00031* 

The  pressure  distribution  along  the  tunnel  was  controlled  by  the 
flew  through  the  slotted  wall  and  the  plenum  ohamber  as  controlled  by  the 
downstream  pressure  control  gate.  In  the  early  exploratory  tests  of  the  pro¬ 
gram  a  slotted  wall  was  used  which  had  a  relatively  high  porosity.  During 
the  tests  it  was  observed  that  flow  passed  through  the  slotted  walls  into  the 
plenum  chamber  in  the  upstream  region  but  in  the  reverse  direction  in  the 
downstream  region.  As  a  result  the  main  stream  velocities  decreased  in  the 
upstream  portion  and  increased  again  downstream.  Correspondingly,  the  pres¬ 
sure  gradient  was  strongly  negative  upstream  and  positive  downstream. 

For  the  tests  documented  in  this  reporta  slotted  wall  of  much  lower 
porosity  was  used.  The  design  of  the  wall  was  in  accordance  with  the  remarks 
made  earlier  on  slotted  walls.  The  porosity  was  made  low  enough  so  that  the 
slot  velocity  approximates  the  core  velocity.  The  porosity  was  also  varied 
in  the  lengthwise  direction  in  proportion  to  the  calculated  variation  in  the 
rate  of  displacement  thickness  increase.  During  a  test  the  pressure  in  the 
plenum  chamber  was  adjusted  until  the  piezometer  tubes  showed  a  minimum  pres¬ 
sure  variation  in  the  tunnel. 

Figure  9  shows  the  readings  for  the  optimum  settings  for  thB  three 
standardised  flews t  2.U  fps,  6.0  fps,  and  15  fps.  The  referenoe  pressure 
used  for  these  plottings  is  the  average  of  all  the  readings  for  any  given 
condition.  The  variation  of  pressure  head  along  the  channel  is  mostly  less 
than  plus  or  minus  1  per  cent  of  the  main  stream  velocity  head.  Without  a 
slotted  wall  the  rise  in  pressure  bead  for  this  length  would  be  about  6  per 
cent  of  the  velocity  head.  The  pressure  head  drop  across  the  slotted  wall 
for  optimum  conditions  was  measured  to  be  0.21*  ft,  0.65  ft,  and  2.Ui  ft  for 
the  2.U,  6.0,  and  15  ffcs  core  velocities.  The  slot  velocities  corresponding 
to  these  heads  are  3.9 ,  6,5,  and  12.2  fps  or  1.62,  1.08,  and  0.81  times  the 
respective  core  velocities. 

Figure  9  presents  pressure  distribution  data  for  the  four  vibration 
conditions  (including  sero  frequency).  The  plenum  pressure  was  adjusted  to 
optimum  for  each  flow  conditionbut  was  kept  the  same  for  all  the  frequencies 
for  each  flew  so  that  any  influence  of  frequency  on  the  pressure  would  be 
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detected.  The  data  Indicate  that  no  measurable  Influence  existed  for  these 
conditions. 

Tbs  most  important  results  to  be  presented  from  this  program  are 
the  boundary  layer  velocity  prof  lies  on  Figs.  10,  11,  and  12.  Figure  10  shows 
the  velocities  taken  at  Stations  1  and  2  for  a  flow  having  a  free  stream 
velocity  of  2.1*  fps.  Figure  11  contains  the  data  for  a  free  stream  velocity 
of  6,0  fps  and  Fig.  12  for  13  fps.  Data  on  the  four  vibrational  conditions 
are  given  for  each  flow  for  Station  2, but  only  zero  frequency  data  are  given 
for  Station  1  upstream.  It  must  be  stated  at  the  start  that  the  results  show 
no  measurable  influence  of  vibration  even  though  the  data  cover  values  of 
c/U^  from  55  to  2.5. 

Conditions  in  the  tunnel  made  it  very  difficult  to  obtain  consis¬ 
tent  data.  The  main  source  of  difficulty  was  a  very  low  frequency,  irregular 
pressure  surge  that  made  manometer  reading  very  dif fic\’ It.  Several  unavoid¬ 
able  factors  seemed  to  contribute  to  this.  These  included  surges  in  the  300 
ft  long  supply  channel,  pressure  pulses  in  the  approximately  100  ft  long 
penstock,  and  air  accumulation  from  air  precipitation  in  the  region  of  sub- 
atmospheric  pressure  Immediately  downstream  of  the  intake  valve.  Deepening 
in  the  nanometer  system  was  ineffective  and  electronic  sensing  devices  that 
were  tried  were  not  accurate  or  consistent  enough. 

The  method  finally  resorted  to  was  to  use  straight  half-inch  glass 
piezometers  and  photograph  the  fluid  levels  fifty  times  for  each  setting  at 
1  sec  intervals  using  a  16  mm  movie  camera.  Each  movie  frame  showed  four 
piezometer  tubes.  The  tube  to  the  left  was  connected  to  the  wall  static  tap 
in  the  tunnel  roof,  located  2  in.  to  the  left  of  the  centerline  at  the  meas¬ 
uring  section;  the  second  from  the  left  was  connected  to  the  boundary  layer 
stagnation  tube;  the  third  was  connected  to  the  main  stream  stagnation  tube 
at  Station  3  and  the  last  was  connected  to  the  wall  static  tap  to  the  right 
of  the  centerline.  The  movie  frames  were  read  on  a  microfilm  viewer  against 
an  arbitrary  scale  attached  to  the  ground  glass.  The  results  were  typed  di¬ 
rectly  onto  punched  tape  for  electronic  computing.  The  computer  calculated 
the  u/U^  ratio  for  each  frame  and  gave  the  mean  value  and  standard  devia¬ 
tion  for  the  fifty  frames  for  each  setting.  Settings  with  suspiciously  high 
standard  deviations  were  rerun  and  the  values  with  the  lowest  standard  devi¬ 
ation  was  used.  With  proper  lighting  the  meniscus  was  very  well  defined  and 
good  results  were  obtained  using  telephoto  lenses  even  for  the  low  flows. 
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The  data  in  Figs.  10,  11 ,  and  12  are  plotted  in  the  most  fundamen¬ 
tal  form  possible  for  several  reasons.  First,  this  plotting  is  sufficient 
for  the  purpose  of  the  report ,  namely  to  demonstrate  the  effect  of  vibration. 
Second,  this  fundamental  font  is  most  convenient  for  others  who  may  wish  to 
replot  it  in  any  way.  Third,  the  dimensionless  plots  require  the  use  of  pa¬ 
rameters  such  as  shear  velocity,  boundary  layer  thickness,  and  distance  from 
the  leading  edge  that  were  not  conducive  to  accurate  measurement.  There  also 
have  been  no  corrections  applied  to  compensate  for  the  presence  of  the  wall, 
the  velocity  gradient,  turbulence  intensity,  and/or  laminar  effects  at  the 
tip  of  the  probe.  All  the  data  shown  were  taken  with  the  smallest  stagnation 
tube,  which  has  a  thickness  of  0.0085  in.  A  correction  of  *  O.OOl*  in.  was 
made  to  the  dial  gage  reading  to  account  for  the  semidiameter. 

The  results  are,  however,  in  fairly  good  agreement  with  that  of 
other  researchers.  In  Fig.  10  the  data  for  y  greater  than  about  0.0U  in. 
can  be  fitted  almost  exactly  to  the  logarithmic  law 

v* 

-V  •  5.85  log  — sL  ♦  5.56  (37) 

V  v 

o 

which  was  suggested  hy  Nikuradse's  measurements  (see  Schlictlng  [6],  page 
U37.)  provided  values  of  V*  of  0.126  and  0.113  fps  are  used  for  Station  1 
and  Station  2  respectively.  For  y  less  than  0.0U  in.  the  data  fall  away 
from  this  curve  and  follow  a  transition  curve  leading  to  the  laminar  sublayer 
(sea  Schlicting  [6],  page  l|05  ). 

Station  1  and  Station  2  are  1.1*6  and  5»80  ft  downstream  from  the 
leading  edge  of  the  tunnel  transition  measured  horizontally.  The  temperature 
of  the  water  during  the  tests  varied  between  1*3°  and  1*5°  F.  The  Reynolds 
number  based  on  distance  from  the  leading  edge,  is  0.25  x  10^  for  Station  1 
and  1.00  x  10^  for  Station  2.  Schlicting  [6],  p.  1*38, presents  values  of  the 
local  drag  coefficient  cf  as  a  function  of  Reynolds  number  for  the  veloc¬ 
ity  profile  in  Eq.  (37) .  The  shear  velocity  easily  follows  and  calculates  to 
be  0.117  and  0.105 fps  for  Stations  1  and  2.  These  values  are  both  very  near¬ 
ly  7  per  cent  lcmer  than  the  ones  obtained  from  the  curve  fitting  process. 
This  is  not  unreasonable  considering  the  state  of  the  art  of  measuring  wall 
shear  at  the  time  these  equations  were  evaluated. 
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Figure  11  shows  the  measured  Telocity  profiles  at  Station  1  and 
Station  2  for  a  free  stream  velocity  of  6.0  fps.  The  velocities  fit  Eq.  (37) 
quite  closely  for  values  of  y  between  about  0.01  and  0.1  in.  when  values 
of  shear  velocity  of  0.285  and  0.251  fps  are  used.  The  Reynolds  numbers  at 
these  stations,  completed  as  above,  are  0.62  x  10^  and  2.U7  x  10^  and  the 
corresponding  values  of  shear  velocity  from  Schlicting  are  0.270  and  0.239 
fps.  The  discrepancy  for  this  flow  is  about  5  per  cent.  For  values  of  y 
greater  than  0.1  in.  the  velocities  are  greater  than  given  by  Eq.  (37) .  In 
this  region  the  velocity  profile  fits  the  l/7th  power  law  much  better. 

Figure  12,  which  contains  the  measured  velocity  profiles  for  a  core 
velocity  of  15  fps,  fits  Eq.  (37)  well  for  y  less  than  0.0U  in.  and  the 
l/7th  power  law  for  y  greater  than  O.OU  in.  The  shear  velocities  obtained 
from  the  curve  fitting  are  0.668  and  0.583  fps  for  Stations  1  and  2.  The 
Reynolds  numbers  are  1.55  x  10^  and  6.20  x  10'b  and  the  corresponding  shear 
velocities  from  Schlicting  are  0.621*  and  0.556  fps. 

The  data  on  Fig.  12 are  the  most  interesting  from  the  point  of  view 
of  the  study  because  the  15  fps  flow  has  the  lowest  wave  celerity  to  core 
velocity  ratio.  The  ratio  has  values  of  2.1*6,  5.70,  and  11.3  for  the  1U.8, 
57,  and  122  cps  vibrations.  These  values  are  within  the  range  where  the  vis¬ 
cous  part  of  the  velocity  disturbance  starts  to  extend  beyond  the  laminar 
sublayer.  The  data,  nevertheless,  show  that  the  mean  properties  of  the  bound¬ 
ary  layer  are  unaffected  by  the  vibration. 

VII.  CONCLUSIONS 

From  the  theory  presented  it  is  possible  to  conclude * 

(1)  A  turbulent  boundary  layer  passing  over  a  smooth  plate 
vibrating  with  low  amplitude  standing  waves  reoeives  a 
harmonic  disturbance  that  can  be  resolved  into  two  addi¬ 
tive  components ,  an  inviscid  motion  and  a  viscous  motion, 
both  of  which  decrease  exponentially  with  distance  from 
the  wall. 

(2)  The  viscous  motion  decreases  very  rapidly  with  distance 
from  the  wall  and  is  largely  contained  within  the  lamin¬ 
ar  sublayer  provided  the  wave  celerity  c  is  greater 
than  3.5  times  the  free  stream  velocity. 
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(3)  The  inviscid  motion  has  very  low  Telocity  gradients  and 
becomes  negligible  at  a  distance  of  0.7A  from  the  wall. 

(U)  A  sheltering  effect  is  possible  if  the  invisoid  motion 
becomes  out  of  phase  with  the  wave  font. 

(5)  The  boundary  condition  of  the  viscous  motion  governs  the 
relationship  of  the  inviscid  motion  to  the  wave  form  and 
is  such  that  there  is  no  sheltering  effect,  provided  that 
c  is  greater  than  3.3  times  the  free  stream  velocity. 

(6)  The  sheltering  effect  if  any  is  &  function  of  the  wave 
oelerityand  the  free  stream  velocity  and  does  not  depend 
directly  on  wave  amplitude ,  wave  length,  or  wave  frequency. 

From  an  inspection  of  the  properties  of  plates  presented  in  the  re¬ 
port  it  is  possible  to  conclude > 

(1)  Stiffness  requirement s  for  any  conceivable  structural  pur¬ 
pose  limit  the  possible  amplitude  of  vibration  to  well 
within  the  applicability  of  the  theory. 

(2)  Stiffness  requirements  also  limit  the  possible  wave  ce¬ 
lerities  to  the  order  of  100  fps  or  greater. 

From  the  experimental  data  it  is  possible  to  conclude  that  standing 
waves  have  no  measurable  influence  on  the  boundary  layer  for  any  conditions 
tested,  even  though  the  vibrating  plate  was  well  over  50  boundary  layer  thick¬ 
nesses  long  and  the  wave  oelerityto  free  stream  velocity  ratio  was  tested  as 
low  as  2.51w 
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Fig.  1  -  Flow  Characteristics  for  a  Stationary  Boundary 
with  a  Moving  Sinusoidal  Wav*  Pattern — No  Main  Flow 
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-  Vibrational  Characteristics  of  a  Flexible  Plate 


Fig.  4  -  Apparatus  for  th»  Vibrating  Plata  Tests 
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Fig.  5  -  Photograph*  of  the  Test  Apparatus 


Fig.  6  -  Boundary  Layor  Pitot  Tuba*  and  Amiably 
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Fig.  7  -  Photograph*  of  the  Boundary  Layer  Pitot  Tube  Alterably 
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Fig.  9  -  Measurements  of  the  Longitudinal  Pressure  Gradient 
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